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By the complex multimode Bogoliubov transformation, we obtain the general forms of squeeze operators 
and squeezed states including squeezed vacuum states, squeezed coherent states, squeezed Fock states 
and squeezed coherent Fock states, for a general multimode boson system. We decompose the squeezed 
operator into disentangling form in normal ordering to simplify the expressions of the squeezed states. 
We also calculate the statistical properties of the SS. Furthermore we prove that if it is non-degenerate, a 
Bogoliubov transformation matrix can be decomposed into three basic matrix, by which we can not only 
check the criterion of the minimum uncertainty state (MUS) found by Milburn Q, but also prove that 
except the special cases, any multimode squeezed state is MUS after the original creation and annihilation 
operators rotate properly. We also discuss some special cases of the three basic matrices. Finally we give 
an analytical results of a two-mode system as an example. 
Oh" PACS numbers: 42.50.Dv, 03.65.-w, 02.30.Tb 
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I. INTRODUCTION 



Squeezed states(SS) are originally important concepts in quantum optics [0. Because of the fascinating property 
that one of their two quadratures has less uncertainty when compared with coherent states, SS have several potential 
, applications, such as in light communication, amplification of signals, and gravitational radiation detection j^]. SS 
C^) also apply to other quantum mechanical systems, such as polaritons and polarons. 

0^ , Early researches focused only one and two-mode SS. Yuen [Q studied one-mode SS generally. Caves and Shumaker 
H discussed two-mode SS in a great detail. When we treat quantum mechanical systems other than quantum optics, 
however, multimode SS inevitably come into sight. Zhang, Feng and Gilmore Q and Lo and Sollie put forward 
• multimode SS, but not in general form in fact. On the other hand, minimum uncertainty states (MUS) is of basic 
,-SZh ' interests in quantum physics. Stoler || found one-mode SS are MUS. And later Milburn Q proved that multimode 
MUS are special multimode SS, which are simply direct products of each one-mode squeezed state. 
' Our present work is to construct the general multimode SS. For convenience, we first give a short review of 
Fock states, coherent states and coherent Fock states in Sec. II. In Sec. Ill we introduces complex Bogoliubov 
transformation (BT) to diagonalize a positive definite bilinear Hamiltonian. Based on the properties of the BT 
matrix, in Sec. IV we construct concisely general multimode squeeze operator (SO), and decompose it into three 
disentangling parts in normal ordering; then we compare our SO with the usual one and two-mode forms. In Sec. 
V we present four kinds of SS, squeezed vacuum states (SVS), squeezed coherent states (SCS), squeezed Fock states 
(SFS), and squeezed coherent Fock states (SCFS), and their reduced expressions; then we derive their covariance 
matrices, average quantum numbers and their fluctuations; we calculate the wave functions of the SS at the base of 
Fock state, coherent state respectively, and also in the coordinate space. In Sec. VI we prove that if the BT matrix is 
non-degenerate, the BT matrix can be decomposed into three basic matrices, each of which is also a BT matrix. Thus 
we can prove that SCS are MUS after the original creation and annihilation operators rotate properly, In addition, we 
analyze some interesting cases of the decomposed matrices, including the case satisfying the criterion for multimode 
MUS given by Milburn jj]]. To demonstrate our theory, in Sec. VII we give an analytical results of a two-mode 
Hamiltonian as an example. Sec. VIII is the conclusion and prospects of the further work. In Appendix A we give 
another approach to obtain the SVS without the knowledge of the SO. 

In this paper index T stands for transpose of a matrix, * for complex conjugate, and f for Hermite conjugate. 
12i (ID means Yli (Ilf) unless stated specifically, and Df is the short notation for an N x N diagonal matrix 
Diag{/i, /2, . . . , /at}, where / can be any symbol except a, since we use D a to denote the Weyl displacement operator. 
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II. FOCK STATES, COHERENT STATES AND COHERENT FOCK STATES 



For a system with Hamiltonian H = Si^»( a ! a i + a i a \)i where dj (aj) is the annihilation (creation) operator for 
the z-th boson, the corresponding eigenstate is the Fock state \n) — \n\) ® IO2) <8> • • • <8> \un)- 

If adding a linear term to the above Hamiltonian, i.e., H = Hq + Kja, + K*aJ, the ground state is a coherent 

state \a) = £> Q |0), where the displacement operator D a = e^i'°'° r0iai ', and ai — — ^f-- The excited state is 
|n;a) = D a \n), which we call the coherent Fock state. 
Fock state wave functions of the above three states are 



(m|n) = n<5 mi ,n,, 
(m\a) =Y\_ e 



1 a 



to,-! 



and 



(m\n; a) = J~| 



\/n%i\ni\ 

Coherent state wave functions of the above three states arc 



min{mi.ni } 

3=0 



*>=n« 



*) =]Je^ Q ^3(l«.l 2 +lftl 2 ) 



and 



n; a 



e /3*a»-i(|a s | 2 + |/3,| 2 )^* _ a *)n< 



— 2~S and Yi 

coordinate representation of the above three states are 



We define two quadrature operators Xj 



2i 



-, and ket \X) = \Xt) <8> \X 2 ) 



7T V^ 7 ^! 



H ni (V2Xi), 



(X\a)=T\(-)h-^-^ 2 + — 

7T 



and 



where H ni is the rij-th Hermite polynomial. 

The covariance matrices of the three states Cov(JT, Y) are 
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D2n+1 Hn 
— Hn D2n+1 



(1) 

(2) 

(3) 

(4) 
(5) 

(6) 

) Xn), then the 
(7) 
(8) 

(9) 



(10) 



The average quantum numbers of the three states (Ni) are 

m, \oti\ 2 , ni + \ai\ 2 . (11) 

The fluctuations of the quantum numbers of the three states (N?) are 

0, M 2 , (2n 4 + l)H 2 . (12) 

It is noticeable that the coherent Fock state has the same covariant matrix as the Fock state. This is because D a 
only displaces the operators by a constant, which has no contribution to the fluctuations. 
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III. BOGOLIUBOV TRANSFORMATION 



Bogoliubov transformation (BT), originally used in superconductivity can diagonalize a definite positive Hamil- 
tonian of a multi-boson system with its annihilation and creation operators in bilinear form. Javannainen |ic| ] 
introduced a general real multimode BT for a weakly interacting Bose-Einstein condensate confined to a spherically 
symmetric harmonic potential. In order to construct complex SS, it is necessary to extend Javannainen's theory to 
complex form. 

Suppose the Hamiltonian of the system is 



= [ a) T a T } 







a 






at 



where both £ and r\ are N x N matrices, and a (a') is a column vector (ai, 0,2 



(13) 



) T . Instead of using symmetric 



Hamiltonian matrix in Ref . [|10| , we arrange the Hamiltonian matrix in Hermitian form for convenience of the latter 
calculations. 

By the hermiticity of the Hamiltonian it is obvious that 



This Hamiltonian can be diagonalized when expressed by a new set of boson operators c,-, c| , i. e. 



'D a 




c 


D a _ 







H = [ C t T c T ] 

The new boson operators relate to the old ones by a BT matrix M as 



c 


= M 


a 




u V 




a 


c t 


at 




V* u* 




_at_ 



where both u, v are N x N matrices. By commutation relations [aj,aj] = [cj,c|] = 5y, [aj,aj] 
[a\, at] = [c|, ct] = 0, we obtain the inverse of BT 

t T 



— 7jt 7J 



KM^K 



(14) 
(15) 

(16) 
and 

(17) 



where K = 



In 
-In 



with /at the TV x iV identity matrix. Therefore M _1 = KM^K, i.e. 
MKM^ = M^KM = K, 



which in block form is 



and 



T T 

uv = vu 



vJu — V T V* 
t T * 



IN 



From the two representations of the Hamiltonian (|l3| ) and (|T^) , we have 



»7* r 



= M f 







M. 



(18) 



(19) 



(20) 



(21) 



Eq. ( ^l| ) multiplied on right by AT Aft J with the consideration of the property of BT matrix M (fl8|) , yields 
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M f = M f 



D n 
-D n 



(22) 



Therefore the problem of diagonalizing the Hamiltonian is simplified to be a familiar algebra problem of solving 
eigenvalues and eigenvectors of a matrix. The eigenvectors are normalized by Eq. (|l8|). 
Define M = D p hM , with the phase matrix 



D, 



D ei * 
D c - 



(23) 



Obviously, just like A77, also satisfies the secular equation (^) and the normalization ( ft8l). T hus for the complex 
BT case, each eigenvector has an undetermined phase. While for the real BT case as in Ref. |10|] , all the eigenvectors 
can be completely determined up to the signs. Further study of such a property and its effects may be very interesting. 
If the Hamiltonian (|l^) is added by linear terms J2i 



H = [ at T a T ] 

it can be diagonalized as 

H = [ c) T c T 
by a new transformation 

= M 
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n* r _ 
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T t 



'D Q " 
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(24) 



(25) 



(26) 



IV. MULTIMODE SQUEEZED OPERATOR 
A. General case 

The unitary SO U should have the following property 



U 



C/ f = M 



So for an arbitrary 2iV x 27V matrix F, we have 



U[a^ T a T ]F 



LTt = [ a t T a T ] M^FM 



Select F satisfying 



M^FM = F, 



(27) 



(28) 



(29) 



Multiply Eq. (J29J) on left by MK, with the use of (|18|), we get KFM = MKF, which means ifF and M arc 
commutative. Thus ifF is a function of M, F = K f(M). 
Consequently the SO can take the form 



[ a^ T a T ] Kf(M) 



,t 



Suppose / = 



/n /12 

/21 /22 



17 = e 

, with all the four blocks being N x N matrices, it can be deduced that 



(30) 
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u 



c/t = 



/n — /12 + /S 

/21 + /Ji /22 - /n 



The form of M in ( fut ) suggests 



12 — /21; /22 — 111- 



And the unitarity of U demands [Kf(M)y = —Kf(M), namely 



r ft _ f t 1 

ft _ft 
j 12 J 22 



" fll ~fl2 
I 21 f22 



These two relations lead to 



So (pit) is reduced to 



fll — — /22, /£ — /l2) fll— hi- 



u 



C/t = e -2/(M) 



Compared with Eq. (p7|), we get 



— InM. 
2 



/(M) 

Therefore the general complex SO of the TV-mode boson system is 



,t 



C/ = e 

It is very interesting that the SO is in the same form in terms of the new set of operators, i. e. 



U = e 



■§ [ C tT C T liClnM 
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(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



This result can be proved directly by using the relation (|l^), but we can easily understand it by the definition of 



the SO (|27|), since it is obvious that U 



W = M 



if we use the relation (ITq) , which is the same form as Eq. 



Similarly, it can be proved that the product of two SO's 



and 



is a new SO as 



Ux = e 



U 2 = e 



-|[ O tT a T ]K\n(M 1 ) 



■i[ a tT a T ]K\n(M 2 ) 



UxU 2 = e 



(39) 



minding the inverse order of the BT matrices Mx and M 2 . 

As best of our knowledge, our SO (|37j) is more general than all previous forms. For example, Lo and Sollie Q| 
constructed a multimode SO 



5 



exp 



Obviously, there are no a%a\, af, and a} 2 terms in the exponent of this form. It results in 

c/t 



u 



cosh |/3| 2i^/3 



sinh|/3| r 



(3* cosh | 





a 




at 



(40) 



where is defined in Ref. 0. Similar result was also reported by Zhang, Feng and Gilmore 0, in which a 2 and a 2 
terms are included. 

In Eq. ( fie) ) the diagonal blocks of the BT matrix are Hermitian, and the two off-diagonal blocks are symmetric, 
which are not the necessary properties of the four blocks defined in expression ( |l6| ) for a general BT matrix. 



B. One-mode and two-mode SO 



It is interesting to compare our SO with the usual one-mode and two-mode ones. For one-mode case, expression 
([}?]) is directly reduced to 

^ (1) -ex P v vi * ; [otT ar] hr - u ; t . (4i) 

1 2^(^)2-1 L v —r~ J L a J I 

Since there is one degree of freedom according to the discussion at the end of Sec. Ill, we can select a real u, and 
thus get 

UW = exp(i(> 2 - \^ 2 ) (42) 

where £ = £M^+H1_ Therefore the usual one-mode SO is indeed a general one. 
We also prove that the usual two-mode SO 

U {2) = expr(e- l¥, aia 2 - e zip a\4) (43) 
can be deduced from expression (|37|) by selecting a special BT with 



cosh r 
cosh r 



e lv sinhr 
e^sinhr 



Therefore the usual two-mode SO is not a general form. 



C. SO in disentangling form 



Suppose the disentangling form of the SO in normal ordering to be 

U = C e~ aiTpaf e -a +T ™-« T °-'a t e <* T ™ (44) 

where Co is a normalization coefficient, p, a and r are all N x N matrices. Furthermore, we can always arrange 
p = p T , cr' = cr T , and r = r T . 
It is straightforward that 
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" 
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ot_ 


e -0 ra = 


2r 


1 




at 



G 



tT T T + 

-a 1 era— a a a' 



„a' era+a a 1 



" e 2CT 
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e- 2ffT _ 







1 2,0 
1 



So 



C/ f = |G)| 2 



ICnl 2 



Comparing this result with (|27j), with the use of (19), we get 



" 1 0" 




' e 2a 




' 1 


2p 
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2t 1 




e~ 2 ° T _ 
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at 


- e 2a 




2e 2a p 






a 




_ 2re 2CT 


4re 2 > + er 2(yT 




_at_ 





I C 1 = 1, P = \ u ly > cr =ilnu, r = ^v*u 1 . 



(45) 



(46) 



When there is no squeezing, U = 1. On the other hand, u = 7jv, w = 0. Thus U = Co according to Eq.'s ( (44|) and 
@.SoC = l. 

This decomposition of SO is very useful in the calculation involving SS, which will be shown in the next section. 



V. FOUR KINDS OF SQUEEZED STATES 

A. Expressions and reductions 

SO is a key to SS, from which we can derive SVS, SCS, SFS and SCFS. Physically, SVS (SFS) is the ground 
(excited) state of the Hamiltonian (E3h, and SCS (SCFS) is the ground (excited) state of that Hamiltonian plus linear 



terms n*cii + jqaj. Respectively, the SVS, SCS, SFS and SCFS are 

\o). = u\o), 

\a) s = e^ a ^- a ' c '\0) s = Ue^ a ' a "- a:a '\0) = U\a), 

\n) s = [] ^%|o) s = cr JJ W|o) = u\m), 



and 



where 



rr. 



Mr 

C/ e E ( a *°I-<«*TTM^|Q) = u\ 



T 



n: a) 



a 


1 




?7 


-l 


K* 


a* 


~~ 2 








K 



(47) 
(48) 

(49) 



(50) 



(51) 



These formal expressions are symmetric for the annihilation and creation operators, and are convenient such as for 
successive different SO's acting on these states, according to the rule (^9|). But when calculating the wave functions 
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which we will do later, such forms are not so convenient. So it is meaningful to use the disentangle formula ( f44| ) to 
simplify the above four expressions, with only creation operators left. 

From ( ff4| ) and (|4^), together with the identity e TrA = |e A |, the SVS can be reduced as 



|0>, 



e 

\u\ 



, g — a era — a a a' ra 

+ T t 

~ a pa |0). 



|0) 



(52) 



It is interesting that we can also get it in a totally different way demonstrated in Appendix A. 

Since the energy of the SVS is TiDq, and the average energy of the original vacuum (0|U|0)is Tr£ =Tr(uDniv + 
v*Dqv t ) =Tr [(vJu + v t v*)Dq\ =Tr[(l + 2v T v*)Dn] >TrDn, we emphasize an important property that the SVS 
are more stable than the original vacuum states. In another word, the original vacuums do not exist in a stationary 
system. This point is obvious physically because it is the SVS other than the original vacuum states that are the 
ground states of the Hamiltonian ( p"3| ) . 

The SCS can be reduced as follows. 



\a) s = Ue aTa '- a ' a U*U\Q) 



[ -a f a T ] m 



= \u*\-*e- a pa 



| u |-5 e -" t7 > + |o) 



1 -2p 
1 



[-at a T] 



|0). 



Using mm and (Ufl), we find 



M 



' 1 


-2p~ 




' u 





1 




* — T 
V u 



So 



t T 



' u 




"at" 


* — T 






V u 




a 



g-^a'a+Q toluol u a' |q\ 



T — T f 



| w |~2 e ^2 a Q + a TQ g Q u a -a* pa 1 |q^ 



|0) 



(53) 



The reduction of the SFS is indirect and a little more complicated. First, 



n 





a 




([«* «*] 


_at_ 


) n * 



J- _„tT„„t 



10) 



\u\-ie-^Yl ^ a + U ^ a ^ |0). 



Since (w*a + u r oi')j and (w*a + it a^j commute for different i and j, 



10) 



~d~^ 

d n * 



a p T (v*a+u r a t ) 



|0) 



10) 



Pi=0 



Pi=0 



where p = [pi,P2, • • ■ ,_Pjv] T - Therefore 
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\n) s = \u\-h- a ' Tpa '' 



11 d p n z 



i-T rp rp X 

P Tp+p U 01 |Q) 



Pi=0 



n 



i ^ 



p p rp+p u a 1 



Based on the above skills, the SCFS can be reduced to 

[ — a T 
{v*a + u- T a))T 



|0>« 



|n;a) s = |it|-*e- atTpa V 



" u 




a 
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10) 



1 +T t T —T t / t T *\ 



Yl 1 c P T rp+p T M Ta '|0) 



Pi=0 



w 2 e 2 



n 



_J_ p T Tp+p T (M~ T a t +2TQ-a') 



10) 



=n 



^_ e V T rp+p T (u~ T a f +2tq-o*) 



Pi=0 



i=0 



(54) 



(55) 



B. Statistical properties 

The covariant matrice Cov(X, Y) of the SVS or SCS is 



1 

4 

Cov(V, V) of the SFS or SCFS is 



(u — vY(u — v) i(u — vy(u + v) 
—i(u + vY (u — v) (u + u)t (u + v) 



(u - v)^ D n+1 (u -v) + (u- v) T D n (u - v)* i(u - u) t D n+ i(u + v) - i(u - v) T D n (u + v)* 
i(u + v) T D n (u - v)* - i(u + v)^D n+1 (u -v) (u + v)^D n+1 (u + v) + (u + v) T D n (u + v)* 



(56) 



(57) 



For the similar reason to the non-squeezing case discussed in Sec. II, where the Fock state and the coherent Fock 
state have the same Cov(X,F), here the SVS and SCS (the SFS and SCFS) have the same Cov(X,Y). 
The average quantum numbers (N?) of the four kinds of state are 



(58) 



and 



(v)D n u)u + (v^D n+1 v)i 



{u^D n u)u + (v^D n+1 v)u + (CC f ) l4 , 



(59) 
(60) 

(61) 



where C = u T a* — v^a. 

The fluctuations of quantum numbers AN? of the four kinds of states are 



„t 



' 1 1 7 



(62) 
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- (u^uiCCF^i - (v^u) u (C*C%, (63) 

(u> D 2n +iv)ii{v^ D 2n +iu)ii + {v)D n+1 u + v 1 D n v) li (v) D n u + v^D n+1 v)u 

~{A + B)^D n(n+1) (A + B)- 2AW n(n+1) B (64) 

and 

(u f D 2 n+iv)u(v i D 2n+ iu)ii + (v)D n+ iu + v^D n v) li {v) D n u + D n+1 v)a 
— (A+ B)^D n{n+1) (A + B) - 2A^D n(n+1) B + (u^D 2n+1 u + «t D 2n+1 v) lt (CC^ u 

- (u^D 2n+1 v) u (CC T ) u - (viD 2n+1 u) u (C*C%, (65) 

where A l j = \uij\ 2 , Bij = \vij\ 2 . 

C. Wave functions 

Start from (|5^-^5|), we can derive the wave functions of the four kinds of SS at the base of Fock state, coherent 
state and coordinate. Some of the results can be compared with the single mode case in Ref. [?]. 
At the Fock state base, the wave functions of the four kinds of states are respectively 

MO). = It*]"* ( II -p=fTW J e~ qTp %=o, (66) 



j ldq j 



(m\a) s = ^-le-laW^ra I Yl -L=^Z- I e-^+^^l^o, (67) 



frrijl dq- 



l -w—w- 1 jrii+mj 

:,„|„>„ = |, t | "J fl^^ e v r P+q u P - q pq] (6g) 

1 11 yJniXmA dp i t dq j ' 



and 



(mm: a) s = \u\ 2 e 2 



a' a+a ra 



-pr 1 d ni+m i 

K i s/nMA dp n *dq m -. 



e P rp+q u p-q pq+p (2ra-a*)+g u a \ p _ q ._ a (69) 



At the coherent state base, the wave functions of the four kinds of states are respectively 

</3|0> a = \u\-h-^"l 3 -^ pl3 \ (70) 

(/3|a) s = V^ 1 ""^"^- 2 ""™^ 1 '^, (71) 



| M |-i e -i /3t ^-^ tp ' 3 



'ini?' w *^ (72) 



and 



•\n;a) s = \ u \~h ^ u ~ la ~hi a '' ( i ~ 2aT Ta + 2 ^ pP*) 



n^=f#i ) e^+f T («- T ^-^+ 2 ™)| Pj=0 . (73) 
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At the coordinate base, the wave functions of the four kinds of states are respectively 



(X\0) s = \u- e -x- 

7T 



(X\a) s = \ u - v \-^£^ e -X T (u-v)- 1 (u+v)X+2X T (u-v)- 1 a-i^ a -i a T (u-vr(u-vr 1 a ^ 



(74) 
(75) 



(X\n) s = \u-v\-h- xT ^~ v ^^ u+ ^ x 



ik 



i d n * 



-ip T («-t ) )*(«-i;)- 1 p+2p T («-«)*X| 

C l?>i=0 



(76) 



and 



(X\n; a) s = \u - v\~? e~ xT ' (u->>r 1 (u+v)x+2X T (u-v^a-^a^ a-± a T («-„)* (u-v)-^ 



— ^p T (iA— u)* (w— 1 p+p T [2(tt— u)* X — (u— u) * (u— v) 1 a — a*] 



e 2 



P;=0- 



(77) 



VI. DECOMPOSITION OF THE BT MATRIX 



A. Decomposition and MUS 



We can constructively prove that a general BT matrix can be expressed by multiplication of three basic matrices. 

From the first equation of (|l9|), Hermitian matrices un* and vv^ can be diagonalized by a unitary matrix S\ 
simultaneously. So we have Sluu^Si — D cosh 2 r and Sluu^Si = D sinh 2 r . The former identity indicates that T\ = 
D i S\u is a unitary matrix, which is to say, u = S f iD cos h T -Ti- Similarly v — SiD s i n \ lr T2- 

cosh r 

Substituting the expressions of u and v into the first equation of (EQ), we can get £> sinh 2 r TiT^' = T1T2 D sinh 2 r . If 
sinh 2 ri values differently for different i, then the unitary matrix TxT^ must be also diagonal, i. e., T2 = D e i V T^. As 
a result, v — SiD e i V sinhr Ti . Then we have 



M 



Si 

si 



D 



e z v sinh r 



d, 



zv sinh r 



Ti 
I? 



(78) 



But we can go a step further to move the factors e lVi of the middle matrix of the r.h.s. of (78) to the other two 
matrices. Since 



-^cosh r D e i<p s i n h r 
D e -i^ s i nnr -^coshr 



D e i V /2 
D -i V /2 



D e -i V /2 
D i V /2 



(79) 



M 



S 

s* 



-^cosh r -^s 



ih r 



ill /' 

cosh r 



T 

T* 



M S M D M T , 



(80) 



where S = SiD eiv ,/2, T = D e - iv ,/ 2 Ti. 

It is easy to check that all of the three factor matrices Ms, Mb and Mj- are BT matrices. Meanwhile, Mg and Mj 
are also unitary and Mjj is Hermitian. 

Consequently the SO is decomposed according to the rule (p9) as 



U = UtU d Us, 



(81) 



where Ut 



■1 -o tT lnTo 
2 c , 



i(a T D r a-a +T _D,.a t ) 



and Us = \S\- 



\ g— a} T In 5a 



Literally, the BT matrix is well known as a symplectic matrix, which is not the ample condition for such a decom- 
position. The other one is that the two diagonal blocks (u and u*) and the two off-diagonal blocks (v and v*) of the 
BT matrix are complex conjugated respectively. 
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There is a problem about the squeezing property of the SVS and SCS. Actually, it is not always squeezed with 



respect to the original quantum a; and a\ according to Table II. But if we define a set of operators 

M T 



b 

&t 



(82) 



then the SVS and the SCS are always squeezed with respect to the new operators, because the covariant matrix with 
respect to the new operators is 



Cov(V, Y) 



1 



D e 2, i 

-i D„- 



(83) 



Unlike Ci, which is a linear combination of both a^'s and aj's, here bi is a linear combination of only a^s, b = Ta. 
Since T is unitary, it means a rotation in the operator space. 

Following Artoni and Birman fllJl , we name the quasi-particles corresponding to the set of bi's and frj's as the mixed 
bosons, and the squeezing as intrinsic squeezing. 

Therefore we conclude that both the SVS and SCS are MUS with respect to the mixed bosons. 



B. Some special cases 



In this subsection we discuss some special cases of BT matrix. 
1. If T = I N , then from Table II, for the SVS and SCS, 



(A.?Q ) — ^[{SDcoshr SDgtnhr') {SD CO shr SD s i n h r )]u — ^D e -2r, 



(AY; 2 ) — -[(S-Dcoshr + SDsinhr)^ (SDcoshr + SAanhr)]™ — -^D e 2r ■ 

So the SVS and SCS are MUS with respect to the original operators, which is just the case discussed by Milburn [Q. 
2. If T = D e ^ /2 , then 

(AX 2 ) = -(cosh2ri — cos <pi sinh 2r{), 



(AY; 2 ) = -(cosh2ri — cos (pi sinh 2r{), 



which is the same as the result of the general one-mode case 

3. If T = 5 f , then u = v) , v = v T 
that there are no a;aj terms, which is just the case discussed by Zhang, Feng and Gilmore [|| and Lo and Sollie @] 



u = u< , v = , i. e. M t = M, and we can get U = e^ al s " DrSi a a ^ T sd t s t a 1 *) ^ j t j g no ticeable 



VII. A TWO-MODE EXAMPLE 



In this section we calculate an example of a two-mode system. 
Supposing a system with the Hamiltonian 



H = 45(aia{ + ajai + a,2a\ + a2 a 2) + 36(a^ + a[ z + a\ + a J 2 z ) 



t2 



from which we get 



\Ai{a\a2 - 


- a\a 2 ) — ?>2i(aia2 — 


a\a\), 




'45 7i ' 




' 36 Wi ' 


-7i 45 


V = 


16i 36 
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According to our theory, the Hamiltonian is diagonalized as 



H = 30(c lC [ + c\cl) + \2{c 2 c\ + 4c2), 

and the BT matrix block 





[7 3i 




2 2i 


u = 


3\/l0 \/l0 
bi 1 

3^2 \/2 


V = 


3^10 \/l0 
2i^2 n 
3 U 



It is obtained that 



InM 



Z^i ^ 2 «M3 

i^i ^ 4 

^2 »M3 iyti! 

i^i 3 fii ifii 



where \i\ 
cos^ 1 ^/ 
So 



7i sec 72 esc 73, ^ 2 
i 



7itan7 2 + lnr, yU 3 = -71 sec 72 cot 73, ^4 

~3~ 

V5' 



71 tan 72 — lnr, with 71 



45 ~^ , 73 = 2 tan" 1 a/375, and r = - 



JJ _ e iMi(aJa 2 +a 1 aJ)+i^3(a 1 a 2 +aJaJ)-i^-(ai+aJ 2 )-i^(a|+aJ 2 ) 

The SO in normal ordering is 



3^/5 -±(7al 2 +5ai 2 +6ia\ah -A - j r (2a 2 +3a 2 +2t75a 1 a 2 ) 



11 



where A = [ a{ 



t „t 



InJ^+Ai(7-3>/5) 



Consequently, the SVS is 



The rotation matrix S = 



\0)s = 

VE-i V5+1 4 
2V3 2V3 

V5+1 4 \/5-l 
2V3 2^3 



ln \/^i?-^(7-3V5) 



01 
a 2 



and fj, 



n -l V222V1T-94 
7+3V1T 

\/222v^-94 



1 ( 7a t 2 +5a t2 +6M t a t 
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T 



10} 



1 + — 

2 ^ ^5 



71' 



7T 



2 + ^5 



(84) 



, and the diagonal block D CO shr = 



Dia g{\/ ! + 375' V 3 ~ 57I } ' jDsinhr = Diag V^ + V i ~ ITf}' Consequently, 

_ gitan _1 (V5— 2)(aja 2 +aia2)g-9'(ai— aj 2 )+^2.(a2— a^Jg— itan -1 3+ 2 v ^ (a^a 2 +aia^) 



VIII. CONCLUSION 



The major achievements in this paper are threefold. First we constructed the general SO for a multimode boson 
system by using the properties of the complex BT matrix. The crux of the derivation is that the 2N x 2N matrix M 
should be considered instead of the isolated N x N blocks u and v. Second, we obtained four kinds of SS, which are 
ground or excited states of a bilinear Hamiltonian or the Hamiltonian plus linear terms, and simplified them by using 
the disentangling form of the SO in normal ordering. Third, we proved that any M is the multiplication of three 
special BT matrices if m* has no degenerate eigenvalues, which leads to an important conclusion that generally SVS 
or SCS are MUS with respect to a set of mixed bosons. 

We will study several physical systems using our multimode SS theory in the near future, which would also help us 
learn more properties of general or special multimode system. 

We would like to thank Dr. Zidan Wang and Dr. Shaolong Wan for useful discussions. 
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APPENDIX A: CONSTRUCT THE SVS BY THE COHERENT STATE 



The annihilation and creation operators act on (a\ yielding respectively 



/ d a n/ i / , t / , * 
{— + —){a\, {a\al = {a\a t 



Defining / = (a\0) s , then 



= <a|c|0) s = {a\ua + v^\0) s = [u{— + |) + va*]f 



The above equation multiplied on left by u 1 yields 



l da* 2 1 



The second equation of (19) indicates that 2p = u l v is a symmetric matrix, satisfying 



d 



jk 



Thus from (A3) we get 



(Al) leads to -^-f = — -rf-/, so /o is just a constant irrelevant with Qj 
The overcompleteness relation of the coherent state is 



n 



doLjda* , 

^— HH = i. 



So 



i°>- = n / 



^l«)(«|0) s = /oI[ 



dajda i ^- a i a - a \ pct * ^ a T a i . qv 

27TI 



Divide a into real and imaginary part, i.e. let a = a; + ii/, then the above equation becomes 



|0). = /oII 



dxidyi -i xT y T ] 



1 + p —ip 
-ip 1 - p 



|0>. 



Using the Gaussian integral formula with a little change from Ref. pa 



■ ■ ■ dx 2n - Y%_, XjA^+Y'Z *jJi = | A |-i e iE^ =1 J*V J : 



together with the two identities 



(A8) can be reduced to 



|0) s - f a 



1 + p —ip 
-ip I- p 



-ip 
1-P 



e 



1 + p -ip 
-ip 1 - p 



1 + p — ip 
-ip 1-p 



1 - p ip 
ip 1 + P 



-1 


■ at " 




iat 



10} = foe- 



t T t 

2 pa |0), 
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/o can be determined from the normalization of 1 0) s = 1 . 



Inserting (|Aq), and using (Al) together with (0|a) = e 2 " a , we get 



,<0|0>. = |/o| 2 (0|e- a p a e~ al ^ |0) 



_ |J | 2 TT / dajda* c -a T a" -a T p* a-a' pa' _ ^ 



Set on — Xi + ij/i, the above equation is transformed into 
l = |/o| 2 n/ (^)exp(-[^ y T ] 
This integral can be calculated according to (|A9|) with J = 0. 



1 + p + p* ip* — ip 
ip* — ip 1 — p — p* 



l = |/o 



1 + p + p* ip* — ip 
ip* — ip 1 — p — p* 



Define unitary matrix V = ^= 

-2 



1 1 

—i i 



then 



v 



1 + p + p* ip* — ip 
ip* — ip 1 — p — p* 



V 



1 2p 

2p* 1 



|l-4p>| 



From ( |46| ) and (19), we have 1 — 4p*p — (u T u*) 1 , so 

/o = e^|«|-i 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



So we deduced the same SVS as in Sec. Ill, except for a free phase 
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